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Abstract. Interpolation techniques play a central role in Astronomy, where one often needs to smooth irregularly 
Ch ' sampled data into a smooth map. In a previous article (Lombardi & Schneider 2001, hereafter Paper I), we have 

^3 ' 



considered a widely used smoothing technique and we have evaluated the expectation value of the smoothed 
map under a number of natural hypotheses. Here we proceed further on this analysis and consider the variance 
of the smoothed map, represented by a two-point correlation function. We show that two main sources of noise 
contribute to the total error budget and we show several interesting properties of these two noise terms. The 
expressions obtained are also specialized to the limiting cases of low and high densities of measurements. A 
number of examples are used to show in practice some of the results obtained. 



> 
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O 1. Introduction 
(N 

Raw astronomical data are very often discrete, in the sense that measurements are performed along a finite number of 
' directions on the sky. In many cases, the discrete data are believed to be single measurements of a smooth underlying 
^ ^ field. In such cases, it is desirable to reconstruct the original field using interpolation techniques. A typical example 
' of the general situation just described is given by weak lensing mass reconstructions in clusters of galaxies. In this 
^ • case, thousands of noisy estimates of the tidal field of the cluster (shear) can be obtained from the observed shapes 
I of background galaxies whose images are distorted by the gravitational field of the cluster. All these measurements 
• can then be combined to produce a smooth map of the cluster shear, which in turn is subsequently converted into a 
^ ^ projected density map of the cluster mass distribution. 

One of the most widely used interpolation techniques in Astronomy is based on a weighted average. More precisely, 
, a positive weight function, describing the relative weight of a datum at the position 6 + (j) on the point 9, is introduced. 
. The weight function is often chosen to be of the form w(|0|), i.e. depends only on the separation |0| of the two points 
' considered. Normally, w is also a decreasing function of \4>\ in order to ensure that the largest contributions to the 
interpolated value at 6 comes from nearby measurements. Then, the data are averaged using a weighted mean with 
the weights given by the function w. More precisely, calling /„ the n-th datum obtained at the position 0„, the smooth 
map is defined as 

E„=l MS - On) 



where N is the total number of objects. In a previous paper (Lombardi & Schneider 2001, hereafter Paper I) we have 



evaluated the expectation value for this expression under the following hypothesis: 
— The measured values {/«} are independent random variables with expectation value 

(In) = fiOn) . (2) 

In other words, the {/„} are unbiased measurements of a field f{9). 
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— The positions {0n} are independent random variables with uniform distribution and density p. In Paper I we 
initially considered a fixed number N of positions inside a field of finite area A; then, we took the continuous 
limit letting N go to infinity with p ~ N/A constant. Equivalently, we considered A'^ to be a Poisson random 
variable with average pA: 

P^W=e-M(^, (3) 
and each location 0„ to be uniformly distributed inside A: 

Po{0n) = ^ . (4) 

In Paper I we have shown that 



(/(0))=y /(0'Kff(0-0')dV. (5) 

Thus, (/) is the convolution of the unknown field / with an effective weight Wcs which, in general, differs from the 
weight function w. We also have shown that Wcs has a "similar" shape as w and converges to w when the object 
density p is large; however for finite p, Wcs is broader than w. 

Here we proceed further with the statistical analysis by obtaining an expression for the two-point correlation func- 
tion (covariance) of this estimator. More precisely, given two points 9a and 6b, we consider the two-point correlation 
function of /, defined as 

Cov(/; Oa, 9b) = {n9A)f{9B)) - (/>A))(/>i3)) (6) 

In our calculations, similarly to Paper I, we assume that /„ are unbiased and mutually independent estimates of f(9n) 
[cf. Eq. (|^)]. We also assume that the {/„} have fixed variance cr^, so that 

( [fn - [/m - f{e,n)\ ) - CjH,,„, . (7) 

The paper is organized as follows. In Sect. H we summarize the results obtained in this paper. In Sect. I we derive 
the general expression for the covariance of the interpolating techniques and we show that two main noise terms 
contribute to the total error. These results are then generalized in Sect. ^ to include the case of weight functions that 
are not strictly positive. A useful expansion at high densities p of the covariance is obtained in Sect. |^. Section 
devoted to the study of several interesting properties of the expressions obtained in the paper. Finally, in Sect, 
consider three simple weight functions and derive (analytically or numerically) the covariance for these cases. Four 
appendixes on more technical topics complete the paper. 



2. Summary 

As mentioned in the introduction, the primary aim of this paper is the evaluation of the covariance (two-point 
correlation function) of the smoothing estimator (|^) under the hypotheses that measurements /„ are unbiased estimates 
of a field f{0) [Eq. (||)] and that location measurements {9n} arc independent, uniformly distributed variables with 
density p. Hence, we do not allow for angular clustering on the positions {On}, and wc also do not include the effects of 
a finite field in our calculations (these effects are expected to play a role on points close to the boundary of the region 
where data are available). Moreover, we suppose that the noise on the measurements {/«} is uncorrelated with the 
signal (i.e., that variance cr^ is constant on the field), and that measurements are uncorrelated to each other. Finally, 
we stress that in the whole paper we assume a non-negative (i.e., positive or vanishing) weight function w(0) > 0. 
Surprisingly, weight functions with arbitrary sign cannot be studied in our framework (see discussion at the end of 
Sect.!). 

The results obtained in this paper can be summarized in the following points. 

1. We evaluate analytically the two-point correlation function of f{9), showing that it is composed of two main terms: 

Cov(/; 9a,9b) = {f{9A)f{9B)) - {f {9 a)) {f {9 b)) ^T, + Tp. (8) 

The term T^. is proportional to and can thus be interpreted as the contribution to the covariance from measure- 
ment errors; the term Tp depends on the signal f{9) and can be interpreted as Poisson noise. These terms can be 
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evaluated using the following set of equations: 

QisA,SB)^ f [e-^^^-^W-^^-^W - 1] d^fl , (9) 
Jn 

Y{sa,sb) = e-iip[pQ{sA,SB)] ■ (10) 

nOO nOO 

C{wa.wb)^p' Asa AsBe''^^^-'^^-Y{sA,SB) ^ p'C[Y]{wa.wb) , (11) 
^0 Jo 

= — / AHwA{e)wB{e)c{wA{e),WB{0)) , (12) 

P Jn 

Tp^- f A^e[j{e)YwA{e)wB{e)c{wA{e),WB{e)) 

P Jn 

+ [ d^Oi [ d^92.f{ei)f{e2)wA{ei)wBie2)\c{wA{ei) + WA{02),WBiei) + WB{e2)) 

Jn Jn L 

-CA{wA{ei))CB{wB{e2))] ■ (13) 

In the last two equations we used the notation wa{6) = w{9a — S), and similarly for wb{6)- moreover the two 
functions Ca and Cb can be obtained from the following limits: 

Ca{wa)^- hm wbC(wa,wb) , Cb{wb) ^ - hm waC(wa,wb) ■ (14) 

p wb^oo p wa^oo 

2. We show that the quantity C{wa, wb) of Eq. (pd]), in the limit of high density p, converges to 

r^t \ ^ I P^S2Q P^Sii P^Sq2 

C{Wa,Wb) — 7 TT T + 7 TT7 7 + 7 v77 v7 + 7 77 7T • (l^) 

[p + WAjKP + WB) [p + WAriP + WB) (p + WArip + WB)'^ (p + Wa){P + Wb^ 

where Sij are the moments of the functions (watWb)'- 

s,,= [ A^e[wA{e)X[wB{e)Y . (i6) 

Jn 

3. We derive a number of properties for the noise terms and the function C{wa, wb)- In particular, we show (1) that 
w A{d)'w B{d)C {w a{(^) ^ uiB{d)) < P^ in every point 0; (2) that the measurement error has as upper bound < cr^; 
(3) that the same error has as lower bound the convolution J WcSA{S)wceB{S) of the two effective weights 
WcSa{0) = wa{0)C{wa{0), oo) and WcSb{0) = wb(^)C'(oo, wb(0)) (cf. [Lombardi fc Schneider 2001| ); (4) that the 
measurement noise converges to Tg. ~ at low densities (p 0) and to T^^ ~ a"' Sn/ p at high densities {p oo). 

3. Evaluation of the covariance 

3.1. Preliminaries 

Before starting the analysis, let us introduce a simpler notation. In the following wc will often drop the arguments 
9a and Ob in CoY{f\0 a,0 b) and other related quantities. Note, in fact, that the problem is completely defined with 
the introduction of the two "shifted" weight functions WAiO) = w{Oa — d) and WB{d) = w{9b — 9). We also call 
/a = /{Sa) and /s = /{Ob) the values of / at the two points of interest 9a and 6b, so that 

= 7T— • (17) 

Hence, Eq. ^ can be rewritten in this notation as 

Cov{f) = {fAfB)~{fA){fB) ■ (18) 

Note that, using this notation, we are not taking advantage of the invariance upon translation of w{0) in Eq. (^); in 
other words, we are not using the fact that wa and wb are basically the same function shifted by 6a — 6b. Actually, 
all calculations can be carried out without using this property; however, we will explicitly point out simplifications 
that can be made using the invariance upon translation. 

We would also like to spend a few words about averages. Note that, as anticipated in Sect. we need to carry 
out two averages, one with respect to {/«} [Eqs. (||) and (^], and one with respect to {6n} [Eqs. (^) and (^)]. Taking 
{6n} to be random variables is often reasonable because in Astronomy one does not have a direct control over the 
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positions where observations are made (this happens because measurements are normahy performed in the direction of 
astronomical objects such as stars and galaxies, and thus at "almost random" directions); it also has the advantage of 
letting us obtain general results, independent of any particular configuration of positions. Note, however, that taking 
{On} to be independent variables is a strong simplification which might produce inaccurate results in some context 



(e.g., in case of a direction dependent density, or in case of clustering; see Lombardi et al. 2001). Finally, since the 
number of observations N is itself a random variable, we need to perform first the average on {/«} and then the one 

on {Ojn}- 

In closing this section, we observe that in this paper, similarly to Paper I, we will almost always consider the 
smoothing problem on the plane, i.e. we will assume that the positions {9n} are vectors of R^. We proceed this way 
because in Astronomy the smoothing process often takes places on small regions of the celestial sphere, and thus on 
sets that can be well approximated with subsets of the plane. However, we stress that all the results stated here can 
be easily applied to smoothing processes that takes places on different sets, such as the real axis M or the space R'^. 

3.2. Analytical solution 

Let us now focus on the first term on the r.h.s. of Eq. (Eq). We have 



(hfB) = ^ I d^e, I d^... . . / d^^.i%^i^^:i^^ . (19) 

A" Jn Jn Jn '^a(0«)J [L™ "'s(^m)J 

Note that the average in the r.h.s. of this equation is only with respect to {/«}• Expanding the numerator in the 
integrand of this equation, we obtain N'^ terms, N of which have n = m and N(N — 1) have n ^ m. We can then 
rewrite Eq. ( p^ ) above as 

{fAfB)=Ti+T2, (20) 

where 

Jn Jn Jn [En ^a(0„)] [E™ ^sl^™)] ^ ' 

T2.^[ d^e. j d^02. . . / d^.. ^f--"^^;;">f"^7;^";^ . (22) 

A" Jn Jn Jn LE„ w^a(^«)J [Em "^sl^™)] 

Despite the apparent differences, these two terms can be simplified in a similar manner. Let us consider first Ti. Using 
Eq. (^, we can evaluate the average (/^) = cr^ + [/(0„)] . Since the positions {On} appear as "dummy variables" in 



Eq. (21]), we can relabel them as follows 
T - ^ f d^f) f d'O ( d^B [/'(^i)+^']^-4(ei)z^B(gi) 

- 7/ ' Vo ^ ■ ■ in ' [E„-A(^„)][Em-B(0m)] ■ ^''^ 

In order to simplify this equation, we use a technique similar to the one adopted in Paper I. More precisely, we split 
the two sums in the denominator of the integrand of Eq. (p3|), taking away the terms ^^^(^i) and wb{Oi). Hence, we 
write 

Ti = - [ d^e^lfiOi) + a^]wAiOi)wB{Oi)C{wA{Oi),WB{Oi)) , (24) 
P Jn 

where C{wa, wb) is a corrective factor given by 

C{WA. wb)^^J d-er ■ ■ / dH^——-^ — - . (25) 

^ Jn Jn [WA + E«=2^^(^")J F-B + Em=2^'-B(^n)J 

The additional factor p = N jA has been introduced to simplify some of the following equations. Note that in the 
definition of C wa and wb are formally taken to be two real variables (instead of two real functions of argument 0\). 
The definition of C above suggests to define two new random variables yA and yB'- 

N 

yx = Y,wx {On) , with X = {A, B} . (26) 
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Note that the sum runs from n = 2 to n = A'^. If we could evaluate the combined probability distribution function 
Py{yA, Vb) for Va and i/b, we would have solved our problem: In fact we could use this probability to write C{wa,wb) 
as follows 

CiwA, wb) = r dvA r dvB , ly^y^;y^\ . (27) 
Jo Jo {WA + yA){WB +yB) 

To obtain the probability distribution Py(jjATyB)i we need to use the combined probability distribution Pw{watWb) 
for wa and wb- This distribution is implicitly defined by saying that the probability that wa^O) be in the range 
[wa,wa + dwA] and wb{0) be in the range [wb,wb + ^wb] is Pu,(wA, wb) dw^ dws. We can evaluate Pw{wa,wb) 
using 

pUwa,wb) ^ J f s{wA~WAie))s{wB-WBie))d^e . (28) 

Turning back to (i/A,yB), we can write a similar expression for py: 

Py{yA, yB) = Ji^ll-l I '^^^2- • ■ / d^eN5{yA - WA2 WAN)5{yB - Wb2 wbn) , (29) 



n Jn 



where for simplicity we have called wxn — wx {On)- Note that inserting this equation into Eq. (|27l ) we recover Eq. (p5[), 
as expected. Actually, for our purposes it is more useful to consider yx to be the sum of N random variables {wxn}- In 
other words, we consider the set of couples {^{wAm WBn)^ ^ made of the two weight functions at the various positions, 
as a set of N independent two-dimensional random variables {watWb) with probability distribution Pw{wa,wb)- 
[Hence, similarly to Eq. (|2^), we consider the weight functions wx to be real variables instead of real functions; the 
independence of the positions On then implies the independence of the couples {wAmWBn)-] Taking this point of view, 
we can rewrite Eq. ( p9| ) as 

Py{yA,yB)= / dwA2 / dwB2 Pw{wA2, Wb2)- ■ ■ / dwAN / dwBN Pw{wAN , Wbn) 

Jo Jo Jo Jo 

xS{yA-WA2 WAN)S{yB Wbn)- (30) 

It is well known in Statistics that the sum of independent random variables with the same probability distribution 



can be better studied using Markov's method (see, e.g., Chandrasekhar 198£; see also Deguchi fc Watson 1987 for 



an application to microlensing studies). This method is based on the use of Fourier transforms for the probability 
distributions Pw and py. However, since we are dealing with non negative quantities (we recall that we assumed 
w(9) > 0), we can replace the Fourier transform with Laplace transform which turns out to be more appropriate in for 
our problem (see Appendix^ for a summary of the properties of Laplace transforms). Hence, we define W{sa, sb) and 
Y{sa,sb) to be the Laplace transforms oi Pi^(wa,wb) and Py{wA,WB), respectively. Note that, since both functions 
Pw and Py have two arguments, we need two arguments for the Laplace transforms as well: 

/"CO />OC 

W{sa,sb) = C[p^]{sa,sb)^ dwA du.Be-™-^^"'>,„(«;A, ws) , (31) 

Jo Jo 

poo pOG 

Y{sA,SB) = C[py]{sA,SB)= dvA dyBe~'^y^-'^y^Py{yA.yB) . (32) 

Jo Jo 



We use now in these expressions the Eq. (28) for p^ and Eq. (|30| ) for py, thus obtaining 

W{sA, sb) = \ [ e-'^AWA{0)-SBWs(e) ^2g ^ (33) 
^ Jq. 

N N -I 

[W(sa,sb)\'''' ■ (34) 



Y{SA, sb) 



AN- 



1 d^e^-- 


■ / d^6'Arexp 


In 


Jq. 



-SA ^ WAn - S-B ^ WBn 
n=2 m=2 



Hence, Py can in principle be obtained from the following scheme. First, we evaluate W{sa,sb) using Eq. (^), then 
we calculate Y{saj sb) from Eq. (^), and finally we back-transform this function to obtain py{yAj yB)- 

Actually, another, more convenient, technique is viable. Following the path of Paper I, we now take the "continuous 
limit" and treat TV as a random variable. As explained in Sect. 0, we can take this limit using two equivalent approaches: 

— We keep the area A fixed and consider iV to be a random variable with Poisson distribution given by Eq. (^. We 
then average over all possible configurations obtained. 

— We take the limit N ^ oo taking the density p = N/A fixed. 
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The equivalence of the two methods can be shown as follows. Let us consider a large area A' D A, and let us suppose 
that the number N' = pA' of objects inside A' is fixed. Since objects are randomly distributed inside A' , the probability 
for each object to fall inside A is just A/A' . Hence N, the number of objects inside A, follows a binomial distribution: 



If we now let N' go to infinity with N' /A' ~ p fixed, the probability distribution for N converges (see, e.g. Eadie et al 
197l|y]o the Poisson distribution in Eq. 

We will follow here the second strategy, i.e. we will take the limit iV, A oo keeping p — N/A constant. In the 
limit A — > cx) the quantity W{satSb) goes to unity and thus is not useful for our purposes. Instead, it is convenient 
to define 

QisA, sb) ^ I [e-^---W-««-- W - 1] d^e = A[W{sA, sb) - l] . (36) 
Jo. 



This definition is sensible because, this way, Q remains finite for A oo. In the continuous limit, Eq. ( |34| ) becomes 

= eQi^A^'^B) _ (37) 



Y{sa,Sb) = lim 



. pQ{sa,sb) 
N 



In order to evaluate C{wa,wb), we rewrite its definition ( |27| ) as 

C{wa,wb) ^ P^ [ dxA [ dxB^^^^^^—^-^—^, (38) 
Jo Jo xaxb 

where xx = yx + wx and 

Cni{xA,XB) = '^{xa - Wa)H(xb - WB)Py{xA - WA, Xb - Wb) ■ (39) 

Here il{xx — wx) are Heaviside functions at the positions wx, i.e. 

H(a:) = r . (40) 

1 otherwise. 



Note that is basically a "shifted" version of py. Looking back at Eq. (|3q), we can interpret the integration present 
in this equation as a very particular case of Laplace transform with vanishing argument. In other words, we can write 

CiwA,WB) = p^C[C^/xAXB]iO,0) ■ (41) 

Thus our problem is solved if we can obtain the Laplace transform of Cw/xaxb evaluated at sa = sb = 0. From the 
properties of Laplace transform [cf. Eq. (D.7)] wc find 

/■OC /"OC 

C[(w{xA,XB)/xAXB]isA,SB) = / ds^/ ds'g Z^{s'a, s'^) , (42) 

JsA J SB 

where Z,^ is the Laplace transform of C^,: 

Z^{sa,sb) = C[U{SA, sb) = e-'^'^^-'^^^YisA, sb) • (43) 
Combining together Eqs. (^), (^2|), and (^ ) we finally obtain 



CiwA, wb) = p' / dsA / dsB e-^^"'^-^«"'«y(sA, sb) = p'C[Y]iwA, wb) ■ (44) 
Jo Jo 

In summary, the set of equations that can be used to evaluate Ti are 

Q{sa,sb)= [ [e-^^'^-^W-'^^'^^W - 1] d^e , (45) 

Y{sa,sb) = c^p[pQ{sa,sb)] , (46) 

C{WA, wb) = p" / Asa / ds^ e-'^^^-'^^^Y{sA, sb) = p" C[Y]{wa,wb) , (47) 
Jo Jo 

Ti^- f d^e[f{e) + a^]wA{e)wB{e)c{wA{e),WB{e)) . (48) 

P Jn 
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These equations solve completely the first part of our problem, the determination of Ti. 

Let us now consider the second term of Eq. (^), namely T2 [see Eq. ^)]. We first evaluate the average in {/„} 
that appears in the numerator of the integrand of Eq. (^), obtaining (/„/,„) = /(0„)/(0„i) [cf. Eq. (0) with n^m\. 
Then we relabel the "dummy" variables {©«} similarly to what has been done for Ti, thus obtaining 



We now split, in the two sums in the denominator, the terms WA{di) + wa{02) and WB{di) + wb{(^2) and define the 
new random variables 

N 

zx = Y.wx {0n) , with X = {A, B} . (50) 

n— 3 

Again, if we know the combined probability distribution pz{zA, zb) of za and zb our problem is solved, since we can 
write [cf. Eqs. (||) and d^)] 

^2^ ^^^'^^ / d'gi [ d^02f{ei)f{e2)wA{9i)wB{92) I dzA [ AzbP.{za.zb) ^ 



^ Jn Jn Jo Jo wa{Oi) + wa{02) + ZA 

1 . ^ 

X — — . (51) 

WB\Ol) + WB{62j + ZB 

Actually, in the continuous limit, zx is indistinguishable from yx (zx differs from yx only on the fact that it is the 
sum of iV — 2 "weights" instead of N — 1; however, N goes to infinity in the continuous limit and thus yx and zx 
converge to the same quantity). Thus we can rewrite Eq. (|l]) as 

T2= I A^Oi [ d'e2f{,ei)f{e2)wA{ei)wB{e2)c{wA{ei) + WA{e2),WB{ei) + WB{e2)) , (52) 

Jn Jn 

where C is still given by Eq. (|47|). 

Finally, in order to evaluate Cov(/), we still need the simple averages (/a) and (/b). These can be obtained 
directly using the technique described in Paper I, where we have shown that the set of equations to be used is 

Qx(sx)= /[e-^-"'-(^)-l] d^^, (53) 
Jn 

Yx{sx) = exp[pgx(sx)] , (54) 

Cx{wx) = P / dsxe-^-"'-Fx(sx) , (55) 



{fx)= / d?ef{e)wx{e)Cx{wx{e)). (56) 



We recall that in Paper I we called the combination WaSx{S) = wx{0)Cx {wx{6)) effective weight [cf. Eq. (||) in the 
introduction]. Alternatively, we can use the quantities Y{sa,sb) and C{'wa,wb) to calculate the correcting factors 
Ca and Cb- From Eqs. ( p3| ) and ( p^ we immediately find 

Oa(sa) = O(sA,0) , Ya{sa)^Y{sa,0) , (57) 

Qb{sb) = Q{0,sb) , Yb{sb)=Y{0,sb) . (58) 



Then, using the properties of Laplace transforms [cf. Eq. (D.IO)], and comparing the definition of 0(10 a, wb) [Eq. (^ 
with the one of Cx{wx) [Eq. (pBh] we find 



Ca{wa)^- lim wbC{wa,wb) , Cb{wb) ^ - lim waC(wa,wb) ■ (59) 

p wb-^oo p wa^oo 

We now have at our disposal the complete set of equations that can be used to determine the covariance of /. 

In closing this subsection we makes a few comments on the translation invariancc for wx (see Sect. 3.1). Since 
wa{6) and wb{0) differ by an angular shift only, the two functions Qa and Qb are the same, so that Ca coincides 
with Cb- Not surprisingly, the two effective weights WcSA and WcSb differ also only by a shift. 
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3.3. Noise contributions 

A simple preliminary analysis of the Eqs. and (152) allows us to recognize two main sources of noise. In fact, a 



term in Eq. (48) is proportional to u , and is clearly related to measurement errors of /, namely 

^2 



2 n 

T„ = — d?ewA{e)wB{e)c{wA{e),WB{e)) . 

P Jn 



(60) 



Other factors entering Cov(/) can be interpreted as Poisson noise. Hence, we call Tpi = Ti — Ta, 7p2 = 72, and 
Tp3 = {fA){fB), so that the total Poisson noise is Tp = Tpi + Tp2 — Tps. Note that the Poisson noise Tp, in contrast 
with the measurement noise Ta, strongly depends on the signal f{9). 

The noise term is quite intuitive and does not require a long explanation. We note here only that this term is 
independent of the field f{d) because we assumed measurements /„ with fixed variance cr^ [see Eq. (0)]. 

The Poisson noise Tp can be better understood with a simple example. Suppose that f{6) is not constant and let 
us focus on a point where this function has a strong gradient. Then, when measuring / in this point, we could obtain 
an excess of signal because of an overdensity of objects in the region where f{6) is large; the opposite happens if we 
have an overdensity of objects in the region where f{9) is small. This noise source, called Poisson noise, vanishes if 
the function f{9) is flat. 

In the rest of this paper we will study the properties of the two-point correlation function. Before proceeding, 
however, we need to consider an important generalization of the results obtained here to the case of vanishing weights. 



4. Vanishing weights 

So far we have implicitly assumed that both wa and wb are always positive. In some cases, however, it might be 
interesting to consider vanishing weight functions (for example, functions with finite support). We need then to 
modify accordingly our equations. 

When using vanishing weights, we might encounter situations where the denominator of Eq. (|l|) vanishes because 
all weight functions w{9 — 9^) vanish as well. In this case, the estimator f{9) cannot be even defined (we encounter 
the ratio 0/0), and any further statistical analysis is meaningless. In practice, when smoothing data using a vanishing 
weight function, one could just ignore the points 9 where the smoothed function f{9) is not defined, i.e. the points 9 
for which ■w{9 — 9n) = for every n. This simple approach leads to smoothed maps with "holes", i.e. defined only on 
subsets of the plane. Hence, if we choose this approach we need to modify accordingly the statistical analysis that we 
carry out in this paper. 

This problem was already encountered in Paper I, where we used the following prescription. When using a finite- 
field weight function, we discard, for every configuration of measurement points {9n}, the points 9 on the plane for 
which the smoothing f{9) is not defined. Then, when taking the average with respect to all possible configurations 
{9n} of f{9), we just exclude these configurations. We stress that, this way, the averages (/(0)) and \f{9')) of the 
smoothing (|^) at two different points 9 and 9' are effectively carried out using different ensembles: In one case we 
exclude the "bad configurations" for 9, in the other case the "bad configurations" for 9' . 

The same prescription is also adopted here to evaluate the covariance of our estimator. Hence, when performing 
the ensemble average to estimate the covariance Cov(/; 9a, 9b), we explicitly exclude configurations where either /a 
or fs cannot be evaluated. This is implemented with a slight change in the definition oi Py, which in turn implies a 
change in Eq. ( ^6| ) for Y. A rigorous generalization of the relevant equations can now be carried out without significant 
difficulties. However, the equations obtained are quite cumbersome and present some technical peculiarities. Hence, we 
prefer to postpone a complete discussion of vanishing weights until Appendix |^ we report here only the main results. 

As mentioned above, the basic problem of having vanishing weights is that in some cases the estimator is not 
defined. Hence, it is convenient to define three probabilities, namely Pa and Pb, the probabilities, respectively, that 
/_4 and /b are not defined, and Pab, the probability that both quantities are not defined. Note that, because of the 
invariance upon translation for w, we have Pa = Pb- These probabilities can be estimated without difficulties. In fact, 
the quantity fx is not defined if and only if there is no object inside the support of wx- Since the number of points 
inside the support of wx follows a Poisson probability, we have Px — exp(— pTrjf), where ttx is the area of the support 
of Wx- Similarly, calling ttaub the area of the union of the supports of wa and wb, we find Pab = exp(— pTTyius)- 



Using Eqs. ([45|) and (46) we can also verify the following relations 



Pab ^ lira YisA,SB) , 1= lim Y{sa,sb)- (61) 



Pa^ hm Y{sa,sb), Pb = hm y(s^,ss), (62) 



sb^oo 



SB^O^ 
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Appendix ^ better clarifies the relationship between the limiting values of Y and the probabilities defined above. In 
the following we will use a simplified notation for limits, and we will write something like Pa = ^(O"*", oo) for the left 
equation in (|6^). 

The only significant modification to the equations obtained above for vanishing weights is an overall factor in 
Eq. (^), which now becomes 

C{wa,wb) = - p , p — C[Y]{wa,wb) ■ (63) 

I - -TA - Pb + PaB 

The factor 1/(1 — Pa — Pb + Pab) is basically a renormalization; more precisely, it is introduced to take into account 
the fact that we are discarding cases where either Ja or Jb are not defined. Note, in fact, that in agreement with the 
inclusion-exclusion principle, (1 — Pa — Pb + Pab) is the probability that the both Ja and Jb are defined. Since the 
combination (1 — Pa — Pb + Pab) enters several equations, we define 

^ (64) 



1-Pa-Pb + P. 



AB 



Equation (|63|) is the most important correction to take into account for vanishing weights. Actually, there are also 
a number of peculiarities to consider when dealing with the probability py and its Laplace transform Y . Fortunately, 



however, these peculiarities have no significant consequence for our purpose and thus we can still safely use Eqs. (45) 
and (|4^). Again, we refer to Appendix ^ for a complete explanation. 

In closing this section, we spend a few words on weight functions with arbitrary sign [i.e., functions w{6) that can 
be positive, vanishing, or positive depending on 0]. As mentioned in Sect. |^, in this case a statistical study of the 
smoothing (|l]) cannot be carried out using our framework. In order to understand why this happens, let us consider 
the weight function 

u;(0) = (l-|0|2)exp(-|0|2) . (65) 

This function is continuous, positive for |0| < 1, and quickly vanishes for large \0\. Let us then consider separately 
the numerator and denominator of Eq. (|l|). The denominator can clearly be positive or negative; more precisely, the 
denominator is positive for points 6 close to at least one of the locations 9m and negative for points 6 which are in 
"voids" (i.e., far away from the locations {On})- Hence, the lines where the denominator vanishes separate the regions 
of high density of locations from the regions of low density. Note that, even for very large average densities p, we still 
expect to find "voids" of arbitrary size (in other words, for every finite density p, there is a non-vanishing probability 
of having no point inside an arbitrarily large region). As a result, there will be always regions where the denominator 
vanishes. The discussion for the numerator is similar but, in this case, we also need to take into account the field 
f{9). Hence, we still expect to have regions where the numerator is positive and regions where it is negative but, 
clearly, these regions will in general be different from the analogous regions for the denominator. As a result, when 
evaluating the ratio between the numerator and the denominator, we will obtain arbitrarily large values close to the 
lines where the denominator vanishes. Note also that these lines will change for different configurations of locations 
{On}- In summary, if the weight function is allowed to be negative, the denominator of Eq. (|l|) is no longer guaranteed 
to be positive, and infinities are expected when performing the ensemble average. 



5. Moments expansion 

In most applications, the density of objects is rather large. Hence, it is interesting to obtain an expansion for C{'WAi wb) 
valid at high densities. 

In Paper I we already obtained an expansion for Ca{wa) (or, equivalently, Cb{wb)) for p oo: 

Ca{wa) ^ — + 1— 73 - -1— 74 + — 1— 75 — ■ (66) 

p + WA (P + wa)-^ (p + wa)^ (p + wa)^ 

In this equation, Sij are the moments of the functions (wa,wb), defined as 

= f d^0[wA{0)y[wB{0)y - (67) 

Jn 



Clearly, in Eq. ( |66[ ) enter only the moments Sto, since the form of wb is not relevant for CAiwA)- Similarly, the 
expression for Cb(wb) contains only the moments Sqj- Note that for weight functions invariant upon translation we 



have Sij — ^ji- 
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Fig. 1. The moment expansion of C{wa,wb) for 1- 
dimensional Gaussian weight functions wa{x) = wb{x) 
centered on and with unit variance. The plot shows the 
various order approximations obtained using the method 
described in Sect. ^ (equations for the orders n = 3 and 
n = 4 are not exphcitly reported in the text; see however 
Tabic B.l in Appendix H). The density used is p = 1. 



Fig. 2. The function C{wa,wb) is monotonically decreas- 
ing with WA and wb, while wawbC{'wa,wb) (scaled in 
this plot) is monotonically increasing. The parameters 
used for this figure are the same as Fig. ^ Note that, since 



Pa 



Ph 



0, we have lim,. 



,0+ wawbC{wa,wb) 



\im^^^Q+ wawbC^wajWb) = in agreement with 
Eqs. (|2|) and (p3|); moreover wawbC{wa,wb) < = I 
as expected from Eq. (p4|). 



A similar expansion can be obtained for C{wa,wb)- Calculations are basically a generalization of what was done 
in Paper I for C(ii;) and can be found in Appendix ^ Here we report only the final result obtained: 

r'r \ ^ I P^S2o P^Sii P^So2 , . 

C [WA, Wb ) ^ 7 TT r + 7 t^tt r + 7 tttt + 7 77 7^ . (68) 

{P + Wa)(P + Wb) {p + WAy{p + WB) (p + WArip + WB)'' (p + Wa)(P + WBr 

We note that using this expansion and Eqs. ( p9|) we can recover the first terms of Eq. (|6^), as expected. Figure |l| left 
shows the results of applying this expansion to a Gaussian weight. For clarity, we have considered in this figure (and in 
others shown below) a 1-dimensional smoothing instead of the 2-dimensional case discussed in the text, and we have 
used X as spatial variable instead of 6. The figure refers to two identical Gaussian weight functions with vanishing 
average and unit variance. A comparison of this figure with Fig. 2 of Paper I shows that the convergence here is much 
slower. Nevertheless, Eq. (|68|) will be very useful to investigate some important limiting cases in the next section. 



6. Properties 

In this section we will study in detail the two noise terms T„ and Tp introduced in Sect. ^.3| , showing their properties 
and considering several limiting cases. The results obtained are of clear interest of themselves; for example, we will 
derive here upper and lower limits for the measurement error that can be used at low and high densities. Moreover, 
this section helps us understand the results obtained so far, and in particular the peculiarities of vanishing weights. 



6.1. Normalization 

A simple normalization property for C{wa, wb) can be derived, similarly to what we have already done for the average 
of / in Paper I. Suppose that f{0) — 1 and that no errors are present on the measurements, so that = 0. In this 
case we will always measure f{9) — 1 [see Eq. (|^)], so that (/a) = (/s) = Ij {IaIb) = 1, and no error is expected on 
/. This result can be also recovered using the analytical expressions obtained so far. Let us first consider the simpler 
case of non- vanishing weights. 

Using Eq. ( ^ ) and (48), we can write the term Ipi in the case f{0) = 1 as 



/■oc poo /■ 

Tpi=p dsA dsBe'"2('*-*''*«M d26iwA(6')u;s(6>)e-''-*'"^(^)-'*«™«(^) . (69) 
Jo Jo Jo. 



10 Jo Jn 

32 



The last integrand in this equation can be rewritten as d^Q/dsAdsB [cf. the definition of Q, Eq. ([45|)]: 

^Q{SA,SE 

DsaOsb 



Tp,=p I dsA I dsBe^^f^-^-^^^g^^ . (70) 



'0 
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Analogously, for Tp2 we obtain [cf. Eq. (^2|)] 
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Tp2 = p / dsA / dsB e- 



pQisA-Ss) 



POO POO 

' dsA dsBef"^^'-^'"^^ 
Jo Jo 



[ d^9i i(;^(0i)e-"'*"'-*(^i)-'*«"'«(^i) 
Jn 

dQ{sA,SB) dQ{sA,SB) 



d^e2 WB{,02)e 



- S AW a(62)- S bW 3(^2) 



dsA dsB 

We can integrate this expression by parts taking pe^^ [dQ/dsB) — [9exp(pQ)/9sB] as differential term: 

..pQisA.ss) Q{sA,SB) 



Tp2 = p 



ds, 



ds/ 



,0 Jo osaOsb 



(71) 



(72) 



We now observe that the last term in Eq. ( ^2| ) is identical to what we founded in Eq. (|70|). Hence, the sum Tpi + Tp2 



Tpi + Tp2 — 



dsA e 



pQ(sa,sb) 



dQ{sA,SB) 

dsA 



p& 



pQ{sa,sb) 



00 -1 00 



SA=0. 



SB=0 



Y{oo, 00) - r(oo, 0+) - r(o+, 00) + y (o+, 0+) = 1 



(73) 



The last equation holds because, for non-vanishing weights, F(0+, 0+) — 1 and all other terms vanishes [cf. Eqs. (|61 
H)]. Hence, as expected, (/4/s) - Tpi + Tp2 - 1 - 

In case of vanishing weights, we can still use Eqs. (|70| ) and (^) with an additional factor i' [due to the extra factor 
in Eq. (|63|)]. The last step in Eq. (73) thus now becomes 



Jpi + rp2 = i/[r(oo,oo) - y(oo,o+) - y(o+,oo) + y(o+,o+)] = i . 

The last equality holds since now Y does not vanishes for large (s^, sb) [see again Eqs. (|6l]-||)]. 



(74) 



6.2. Scaling 

Similarly to what was already shown in Paper I, for all expressions encountered so far some scaling invariance properties 
hold. 

First, we note that, although we have assumed that the weight functions wa and wb are normalized to unity, all 
results are clearly independent of their actual normalization. Hence, a trivial scaling property holds: All results (and 
in particular the final expression for Cov(/)) are left unchanged by the transformation w(0) 1— > kw{9) or, equivalently, 

WA{e) >^ kwA(e) , WB{e) ^ kwB{e) . (75) 

A more interesting scaling property is the following. Consider the transformation 

w{e) ^ k^w{ke) , p^k'^p, (76) 

where both factors k^ must be changed according to the dimension of the 6 vector space. If we apply this transformation, 
then the expression for Cov(/) is transformed according to 

Cov(/; Oa, Ob) ^ Cov(/; fc0^, fc^s) . (77) 

This invariance suggests that the shape of Cov(/) is controlled by the expected number of objects for which the two 
weight functions are significantly different from zero. Hence, similarly to what done in Paper I, we define the two 
weight areas A a and Ab ss. 



A 



X 



[wx{e)Yd^6 



S20 if^-A, (73) 
S^i iiX = B. 



For weight functions invariant upon translation we have Aa = Ab • We call Mx = pAx the weight number of objects 
(again, J^a = J^b because of the invariance upon translation). Note that this quantity is left unchanged by the scaling 
([76|). Similar definitions hold for the effective weight Wcsx{0) = wx{d)Cx {wx{6)) and the effective number of objects 
J^cSX = pAcsx ■ 
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6.3. Behavior of C 

In order to better understand the properties of C, it is useful to briefly consider its behavior as a function of the 
weights WA and wb- 

We observe that, since Y{sa,sb) > for every {sa,sb) [see Eq. (^)], C{wa,wb) decreases if either wa or wb 
increase. In order to study the behavior of the quantity wawbC{wa,wb) that enters the noise term Ti, we consider 
the quantity waC{wa, wb). 



waC{wa,wb) = vp I dss 



Ha [wn- Ta f dY(sa,SB) \ 



(79) 



This equation can be shown by integrating by parts the integral over sa- The partial derivative required in Eq. (79) 
can be evaluated from Eq. (^): 

dYjsA, sb) _ ^ dQjsA, sb) ^ms^.ss) ^ (80) 
dsA dsA 

Since this derivative is negative, we can deduce that the integral over sa in Eq. ( [79| ) increases with wa, and thus 
'waC{wa,wb) also increases as wa increases. Similarly, it can be shown that wbC{wa,wb) increases as wb increases. 
In summary, the quantity wawbC{wa,wb) behaves as wawb, in the sense that its partial derivatives have the same 
sign as the partial derivatives of wawb (see Fig. Also, since C{wa,'wb) decreases if either wa or wb increase, we 
can deduce that wawbC{watWb) is "broader" than wawb- 

Since C(wA{d),WB{d)) is positive, the function WA{d)wB{d)C(wA{d),WB{d)) shares the same support as 
WA{S)wB{d). It is also interesting to study the lim its of wau!bC{watWb) at high and low values for wa and wb- 
From the properties of Laplace transform [sec Eqs. ( D.10| )], we have 

lim wawbC{wa,wb) ^ vp^ lim Y{sa,sb) ^ vp^Pab ■, (81) 

0-1- S A — OO 

- T SB^OO 

where Eq. ( |6l| ) has been used in the second equality. Hence, the quantity wawbC{wa, wb) goes to zero only if Pab ~ 0. 
In other cases, we expect a discontinuity at wa = wb =0. Similarly, using Eqs. (|6l|~|6^) we find 

lim wawbC{wa,wb) — v(? lim Y^sa^sb) — v(?Pa , (82) 



-0 



+ 



SB — *00 



lim wawbC{wa,wb) — i^P ^^^^^Y{sa, sb) ^ i^p Pb , (83) 



+0 



+ 



lim wawbC{wa,wb) — i^p hm Y{sa,sb)~i^P . (84) 

WA ^oo 



SA^O^ 

^0+ 



Since wawbC(wa, wb) increases with both wa and wb, the last equation above puts a superior limit for this quantity: 
wawbC{wa,wb) < vp^ . (85) 



6.4. Large separations 

Suppose that the two points Oa and 0b are far away from each other, so that WA{d)wB{S) is very close to zero 
everywhere. In this situation we can greatly simplify our equations. 

If Oa is far away from Ob, then wa{0) and wb{0) are never significantly different from zero at the same position 
0. In this case, the integral in the definition of Q{sa, sb) [see Eq. (|4^)] can be spfit into two integrals that corresponds 
to Qa and Qb [Eq. (|5|)]: 

Q{sa,sb) ^ QAisA) + Qb{sb) , Y{sa,sb) ^Ya{sa)Yb{sb) , C{wa,wb) ^ Ca{wa)Cb{wb) ■ (86) 

Hence, if the two weight functions wa and do not have significant overlap, the function C{waiWb) reduces to the 
product of the two correcting functions Ca and Cb. 

In general, it can be shown that C{wa,wb) > Ca{wa)Cb{wb)- In fact, we have 

/•OO /"OO 

C{wa,wb)-Ca{wa)Cb{wb)^ p"^ dsA dsB e-'^"'^-'^'^^ [eP'^^'-^^"^'' - eP'^^^'-^'f+P'^^^'^''] . (87) 

Jo Jo 
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We now observe that 

Q{sA, SB) - QaIsa) - Qb{sb) = / [e-^-*™^(^) - 1] [e-'B^-(^) - l] > . (88) 

Hence, Q{sa, s_b) > Qyi(sA) + Qb{sb) and the difference between the two terms of this inequahty is an indication of 
overlap between the two weight functions wa and wb - Since the exponential function is monotonic, we find Y{sa, sb) > 
yA{sA)YBisB) and thus 

C{wA, wb) > Ca{wa)Cb{wb) ■ (89) 



6.5. Upper and lower limits for T„ 



The normalization property shown in Sect. 3.1 can also be used to obtain an upper limit for T„. We observe, in 
fact, that Tg. is indistingu isha ble from cr^Tpi for a constant function f{9) = 1. This case, however, has already been 



considered above in Sect. 6A: There we have shown that Tpi + Tp2 = 1. Since Tp2 > 0, we find the relation To- < cr^. 

The property just obtained has a simple interpretation. As shown by Eq. (|6^), is proportional to 1/p and thus 
we would expect that this quantity is unbounded superiorly. In reality, even when we are dealing with a very small 
density of objects, the estimator (^ "forces" us to use at least one object. This point has already been discussed in 
Paper I, where we showed that the number of effective objects, AfcS, is always larger than unity. The upper limit found 
for To- can be interpreted using the same argument. Note that this result also holds for weight functions with finite 
support. 

A lower limit for T^, instead, can be obtained from the inequality (^): 

> — / WAie)wB{e)CA{wA{e))CB{wBie)) d^e^— w,sA{e)w,sB{e)d^e . (90) 

P Jn P Jn 

Hence, the error is larger than a convolution of the two effective weight functions. In case of finite-field weight 
functions, the limit just obtained must be corrected with a factor i/. The argument to derive Eq. ( po|) is then slightly 
more complicated because of the presence of the Px probabilities. However, using the relation PaPb < PAB^ we can 
recover Eq. (|9^) with the aforementioned corrective factor. 

6.6. Limit of low and high densities 

In the limit p — )■ we can obtain simple expressions for the noise terms. If p vanishes, we have Y{sat sb) = 1 [cf. 
Eq. d^)] and thus 

2 

C{wa,wb)^^^ . (91) 
wawb 



In this equation we have assumed wawb > 0. Note that we have reached here the superior limit indicated by Eq. (pSf). 
In the same limit, p 0, Px — 1 — ttxP and v ~ ^/pT^AnB, where tt^hs = t^a + t^b — t^Aub is the area of the 
intersection of the supports of wa and wb ■ Hence we find 

CiwA: wb) =^ . (92) 

An B wawb 

Analogously, in the same limit, we have found in Paper I 

Cx{wx) 0^ , (93) 

TTXWX 

where wx > has been assumed. We can then proceed to evaluate the various terms. For Ta we obtain the expression 
T^^ — I — ^ d^e - . (94) 

P Jr^AHB ^AnB 

Note that the integral has been evaluated only on the subset of the plane where wawb > 0; the case where this product 
vanishes, in fact, need not to be considered because the quantity wawbC{wa,wb) vanishes as well. Exactly the same 
result holds for weight functions with infinite support. Hence, when p — > we reach the superior limit discussed in 



Sect. |6J for T^. 

Equation (|9^) can be better appreciated with the following argument. As the density p approaches zero, the 
probability of having two objects on ttaub vanishes. Because of the prescription regarding vanishing weights (cf. 
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beginning of Sect. the ensemble average in our limit is performed with one and only one object in ttahb- Since we 
have only one object, this is basically used with unit weight in the average (^^, and thus the measurement noise is 
just given by T„ = . 

Let us now consider the limit at low densities of the Poisson noise, which, we recall, has been split into three terms, 



(95) 



Tpi, rp2, and Tpg (see Sect. 3.3). Inserting Eq. (63) into Eq. (04), we find for Tpi 



Tpi 



in 



where (/^)^^ denotes the simple average of on the set ttahb- Hence, Tpi converges to the average of on the 
intersection of the supports of wa and wb- Again, we can explain this result using an argument similar to the one 
used for Eq. ( p^). Regarding Tp2 = T2, we observe that this term is of first order in p because C{wa, wb) is of first 
order [cf. Eqs. (|92|) and (|5^)]. We can then safely ignore this term in our limit p ^ 0. Finally, as shown in Paper I, at 
low densities the expectation value for fx is a simple average of / on the support of wx, i-e- {fx) — {f)nx- Hence, 
Tp3 = (/)7rA(/)7rB and the Poisson noise in the limit of small densities is given by 

2\ 



T, 



{D^aU)^ 



(96) 



In case of a constant function f{6), this expression vanishes as expected. Surprisingly, in general, we cannot say that 
Tp > 0. Rather, if wa 7^ wb, and if in particular the two weight functions have different supports, we might have 
a negative Tp. Suppose, for example, that / vanishes on the intersection of the two supports Trains i but is otherwise 
positive. In this case, the first term in the r.h.s. of Eq. ( p6| ) vanishes, while the second term contributes with a negative 
sign, and thus Tp < 0. On the other hand, if wa = wb then Tp has to be non-negative. 

We now consider the opposite limiting case, namely high density. In this case, it is useful to use the moment 
expansion (|68|). Since T^ and Tpi have an overall factor 1/ p in its definition [cf. Eq. (|60|)], we can simply take the 0-th 
order for C{wa,wb), thus obtaining 



2 /> 

T, ~ — / WA{e)wB{e)d^B^ 
p J 



T 



PI 



WA{e)wBie)f^ie)<f9 



(97) 



For Tp2 and Tp3, instead, we need to use a first order expansion in l/p for C{wa,wb)- This can be done by using 
the first terms in series (|66|), and by expanding all fractions in terms of powers of Inserting the result into the 
definitions of Tp2 and Tp3 we obtain 



Wa{Oi) Wa{02) Wb{0i) ■Wb{02) , S20 + Sii+Si, 



wa{Oi 



Tp2~ /d^^i f d^e2f{0i)f{02)wA{ei)wB{e2) 
Tp3~ f d^9i f d^e2f{e,)f{e2)wA{ei)wB{e2) 

Jn Jn 

Note that we have dropped, in these equations, terms of order higher than 1/p. The difference Tp2 — Tp3 is 

Tp2-Tp3~- f d^Oi j d^e2f{e^)f{e2)wA{ei)wB{e2)[Sii-WA{e2)-WB{9i)\ . 

P Jn Jn 

Using Eqs. (|100|) and (|9l), we can verify that Tp vanishes if / is constant, as expected: 



02 



P P 

Wb{02) >5'20 + Sq2 
P P 



(98) 
(99) 

(100) 



Tpi + Tp2 — Tp3 



P 



1 

P Jn 



dOi / d02 [SllWAiei)wB{e2) - Wa{9i)wa{92)wb{92) - WA{ei)wB{e2)wB{9i)\ 



p P 



(101) 



where the normalization of w has been used. Also, it is apparent that all noise sources, including Poisson noise, are 
proportional to 1/p at high densities. 

In order to further investigate the properties of Poisson noise at high densities, we write it in a more compact form. 
Let us define the average of a function g{6) weighted with q{6) as 



{9), 



g{e)q{e)dH 



q{e)d^ 



Using this definition we can rearrange Eqs. ( |97| ) and (100) in the form 



T 



P 



if). 



if) if) 1 



(/)».)((/>. 



(102) 



(103) 
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This expression suggests that the Poisson noise is actually made of two different terms, ~ (Dwawb (Dwawb 

and {{f)wAWB — {f)wA){{f)wAWB — {I)wb)- The first term is proportional to the difference between two averages of 
and /; both averages are performed using wawb as weight. Hence, this term is controlled by the "internal scatter" 
of / on points where both weight functions are significantly different from zero; it is always positive. The second term 
is made of averages / using different weight functions. It can be either positive or negative if wa 7^ wb- Actually, as 
already seen in the limiting case p — s- 0, the overall Poisson noise does not need to be positive, and anti-correlation 
can be present in some cases. 



6. 7. Limit of high and low frequencies 

The strong dependence of the Poisson noise on the function f{0) makes an analytical estimate of this noise contribution 
extremely difficult in the general case. However, it is still possible to study the behavior of Tp in two important limiting 
cases, that we now describe. 

Suppose that the function f(6) does not change significantly on the scale length of the weight functions WA{d) and 
wb{0) (or, in other words, that the power spectrum of / has a peak at significantly lower frequencies than the power 
spectra of wa and wb)- In this case, we can take the function / as a constant in the integrals of Eq. (fsh, and apply 



the results of Sect. 3.1. Hence, in the limit of low frequencies, the Poisson noise vanishes. 

Suppose now, instead, that the function f{0) does not have any general trend on the scale length of the weight 
functions, but that instead changes at significantly smaller scales (again, this behavior is better described in terms of 
power spectra: We require here that the power spectrum of / has a peak at high frequencies, while it vanishes for the 
frequencies where the power spectra of wa and wb are significantly different from zero). In this case, we can assume 
that integrals such as 

f{e)wx {e)A^e~Q {A, b} (io4) 

n 

vanish approximately, because the average of / on large scales vanishes (remember that we are assuming that / has 
no general trend on large scales). Similarly, the integrals that appear in Tp2 and Tp^ vanish as well. In this case, then, 
the only contribution to the Poisson noise arises from Tpi. This can be easily evaluated 

Tp ~ Tpi ~ / d^ewAi9)wBi9)C{wAie),WB{9)) , (105) 

where we have denoted with the average of |/p on large scales. Hence we finally obtain 

Tp ~ ^^T^ . (106) 

The results discussed in this section can also be numerically verified in simple cases. Figure ||, for example, shows the 
Poisson noise expected in the measurement of a periodic field when using two Gaussian weight functions (see Sect. 
for details). From this figure, we see that the Poisson noise increases with the frequency of the field /, and quickly 
attains a maximum value at high frequencies. Moreover, the same figure shows that, in agreement with Eq. ( |106| ), the 
Poisson noise at the maximum is simply related to the measurement noise T„ (cf. Fig. |^ for p ~ 2). 



7. Examples 

Similarly to what has been done in Paper I, in this section we consider three typical weight functions, namely a top-hat, 
a Gaussian, and a parabolic weight. For simplicity, we will consider 1-dimensional cases only; this will have also some 
advantages when representing the results obtained with figures. Hence, we will use x instead of 6 as spatial variable. 



7.1. Top-hat 

The simplest weight that we can consider is a top-hat function, defined as 

'^i^'*."^- (107) 

otherwise . 



Since w is either 1 or 0, we just need to consider C(l, 1) to evaluate To-. Regarding the Poisson noise, from Eq. (52) 
we deduce that C(l, 2), C(2, 1), and C(2, 2) are also required. 
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Fig. 3. The value of C(l,l)/p for top-hat weights as a Fig. 4. The noise term for two top-hat weights as a 
function of the density p. Both weight functions wa and function of the point separation 5 — \xa ~ xb\ for two 
wb are top-hats [see Eq. (107)] centered on zero. Using densities, p — 2 and p — 5. The plot also shows the quan- 
Eq. (108 ), we can use this graph to obtain Tg- as a function tity Sn/ p, which at high densities approximates T^- (since 
of the density and the point separation xa — xb- then C ~ 1). Note that 5*11 for a top-hat function is just 

given by = (2 - S)/A. 



Figure || shows C(l,l) and C(l,l)/p as functions of the density p for two identical top-hat weight functions 
centered on the origin. From this plot we can recognize some of the limiting cases studied above. In particular, the fact 
that C(l, l)/p goes to unity at low densities is related to Eq. (92); similarly, the limit of C(l, 1) at high densities is 
consistent with Eq. (|68|). The same figure shows also the moments expansion of C(l, 1) up to forth order. As expected, 
the expansion completely fails at low densities, while is quite accurate for p > 5. 

Curves in Fig. || have been calculated using the standard approach described by Eqs. (^5|), (^ ) and (|63|). Actually, 
in the simple case of top-hat weight functions, we can evaluate C(l, 1) using a more direct statistical argument. We 
start by observing that in our case, for xa — xb, we have 

= (72^(1, l)/p. (108) 

On the other hand, a top-hat weight function is basically acting by taking simple averages for all objects that fall 
inside its support. This suggests that, for xa = xb, we can evaluate its measurement noise as 



N=l 



pjN) 
N 



(109) 



where p{N) is the probability of having N objects inside the support. This probability is basically a Poisson probability 
distribution with average p. However, since we are adopting the prescription of "avoiding" weight functions without 
objects in their support, we must explicitly discard the case N = and consequently renormalize the probability. In 
summary, we have 



piN) 



e-Pp^ 

m 



[1 



This expression combined with Eq. (10£) allows us to evaluate C(l, 1) — pT^/a"^ 
e-p pN+i 

N=l 



We can directly verify this result using Eqs. ([451), (46) and (|f 



In fact, for the top-hat function we find 



Q{sA, sb) = 
Y{sa,sb) = 

C(l,l) = 



1] 



k=0 



-k(sA+SB) pk 

fc! 



(110) 

(111) 

(112) 
(113) 



p 

1 - C~P 

e-p 



dsA 

"'0 



dsBY{sA,SB)e~ 



■sa-sb _ 



oo 



p'^e-P 



1-c-P ^ 

k=0 



dsA / dsB- 

-'0 



^{k+l){sA+SB)pk 

kl 



„fe+2 



1-c-P ^fc! (fc4-l)2 



(114) 
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Fig. 5. Numerical calculations for l-dimcnsional Gaussian 
weight functions wa — wb centered on and with 
unit variance. The various curves shows the function 
waWbC{wa,wb) for different densities p. Note that, as 
expected, C{wa, wb) approaches unity for large densities. 
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Fig. 6. Same as Fig. |5|, but for two Gaussian weight func- 
tions centered on and 1 and with unit variance. 
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T„la'\ p = 2 






Sii/p, p = 5 






T,la^ p = 5 _ 


1 








Fig. 7. The noise term for two Gaussian weights (of 
unit variance) as function of their separation. Similarly to 
Fig.lJ, the plot also shows the high-density approximations 
Sii/p. Note that in this case 6*11 is also a Gaussian (with 
double variance). 



Fig. 8. The Poisson noise Tp for two Gaussian weights (of 
unit variance) for a periodic function of the form f{x) ~ 
sin kx as a function of the weight separation d ~ \xa — xb\, 
for a density p — 2. Note that, as expected, the Poisson 
noise increases with k and approaches the limit discussed 
in Sect. 6.7 for high frequencies. More precisely, since for a 
sine function we have (sin^ x) — 1/2, Eq. (106) gives Tp ~ 
T(j/ (2cr^) (this can indeed be verified by a comparison with 
Fig. Note also that, while Tp is strictly positive for 
(5 = 0, it can became negative (see curve for k — 0.5) at 
larger separations. 



Finally, with a change of the dummy variable fc 1-^ n — 1 we recover Eq. (111). 

The other terms needed for the Poisson noise can be evaluated using a calculation similar to the one performed in 



Eq. (114). Actually, it can be shown that for any positive integers wa and we have 



C{wa,wb) 



e-P ^ fc! ik + WA)ik + WB) 



(115) 



Figure ^ shows the expected measurement noise Tcr as a function of the point separation S = \xa — xb\- Note that, 
for densities of order p = 5 or larger, a good approximation is obtained by just taking C(wa, wb) — 1 [cf. the moments 
expansion (|6^)], so that T^. ~ Su/p; we also observe that, for a top-hat weight function, Su is a linear function. 
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Fig. 9. Numerical calculations for 1-dimcnsional parabolic 
weight functions wa = wb centered on and with 
unit variance. The various curves shows the function 
waWbC{wa,wb) for different densities p. 



Fig. 10. The noise term TV for two parabolic weights as a 
function of their separation (see also Figs. ^ and 



7.2. Gaussian 

Frequently, a Gaussian weight function of the form 

w{x) = ^e--V2 (116) 

is used. Although it is not possible to carry out analytical calculations and obtain C^wajWb), numerical integrations 
do not pose any problem. Figure H shows, for different densities, the function wawbC{wa, wb) for two identical weights 
Wa = Wb centered in zero; Fig. shows the same quantity when one of the weight function is centered at unity. Note 
that, in this last figure, the largest covariance is at a; = 0.5, as expected. 

Figure |^ shows the expected measurement noise as a function of the weight separation. Similarly to the top-hat 
weight, an approximation valid for high density is To- = Sn/ p. Figure ^ shows, instead, the Poisson noise Tp expected 
for a field / of the form f{x) = sin/cx, for different values of k . Note that the noise, as expected, increases with k, 
and quickly reaches the "saturation" value discussed in Sect. |6.7| . Note also that the noise is, at lowest lowest density, 
negative for 5 ~ 2.5. 



7.3. Parabolic weight 

Finally, we study of a parabolic weight function of the form 

, , (3x^/4 if bl < 1 , , , 

w{x) = { ' ' ' ' (117) 
I otherwise . 

This function illustrates well some of the peculiarities of finite support weight functions. Figure ^ shows the results of 
numerical integrations for wawbC'{wa, wb) &t different densities p. A first interesting point to note is the discontinuity 
observed at a; = 1, which is in agreement with Eq. (pl[). Moreover, as expected from Eq. (p2[), the function plotted 
clearly approaches a constant at low densities p. Finally, the measurement noise T„ is plotted in Fig. M. 



8. Conclusions 

In this article we have studied in detail the covariance of a widely used smoothing technique. The main results obtained 
are summarized in the following items. 

1. The covariance is composed of two main terms, and Tp, representing measurement errors and Poisson noise, 
respectively; the latter one depends on the field / on which the smoothing is performed. 

2. Expressions to compute To- and Tp have been provided. In particular, it has been shown that both terms can be 
obtained in term of a kernel C'{wa, wb), which in turn can be evaluated from the weight function w{6). 

3. We have obtained an expansion of the kernel C{watWb) valid at high densities p. 

4. We have shown that T^ has an upper limit, given by ct^, and a lower limit, provided by Eq. (pO[). 

5. We have evaluated the form of the noise contributions in the limiting cases of high and low densities. 
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6. We have considered three typical cases of weight functions and we have evaluated C{wa, wb) for them. 

Finally, we note that although the smoothing technique considered in this paper is by far the most widely used 
in Astronomy, alternative methods are available. A statistical characterization of these methods, using a completely 
different approach, will be presented in a future paper (Lombardi & Schneider, in preparation). 

Acknowledgements. This work was partially supported by a grant from the Deutsche Forschungsgemeinschaft, and the TMR 
Network "Gravitational Lensing: New constraints on Cosmology and the Distribution of Dark Matter." 



Appendix A: Vanishing weights 



In Sect. 3.2 we have obtained the solution of the covariance problem under the hypothesis that the weight function w{0) is 
strictly positive. In this appendix we will generalize the results obtained there to non-negative weight functions (see also Sect. 

If wa is allowed to vanish, then we might have a finite probability that j/a vanishes, i.e. a finite probability that no point On 
is inside the support of wa- A finite probability in a probability distribution function appears as a Dirac's delta distribution. 
Since this point is quite important for our discussion, let us make a simple example. Suppose that ^ is a real random variable 
with the following characteristics: 

— ^ has probability 1/3 to vanish. 

— ^ has probability 2/3 to be in the range (0, co); in this range ^ has an exponential distribution. 
Then we can write the probability distribution function for ^ as 

pdo-\m + \m)^M-i), (A.i) 

where H is the Heaviside function [see Eq. (|4o|)]. In other words, the probability distribution for ^ includes the contribution 
from a Dirac's delta distribution centered on ^ = 0. If is known, the probability that ^ is exactly zero (1/3 in this example) 
can be obtained using 

P(C = 0)= f P«(C')dC'= lini /%e(Ode'. (A.2) 

Jo- «-o+ Jo 

Let us now turn to our problem. As mentioned above, for vanishing weights we expect that yA might vanish, i.e. its probability 
might include the contribution from a delta distribution centered on yA = 0; similarly, if wb is allowed to vanish, the probability 
distribution for ys might include a delta centered in ys = 0. For a given ys, the probability PAiys) that yA vanishes is given 
by 

rVA 

PA{yB)= lim / Py{y'A,yB)Ay'A= lim £a [ph(-, ys)] (sa) , (A.3) 

MA-O+Jo- SA^oo 

where the properties of Laplace transform have been used in the last equality (see App endix O) . A similar equation holds for 



the probability that ys vanishes, PsiyA)- Note that the Laplace transform in Eg. (jV.3) is performed only with respect to the 
first variable. The joint probability Pab that both yA and ys vanish is [cf. Eq. (|6l|)] 

rVA rVB 

Pab= 1™ , / / dys py(yA, ys) = ^lim^/:bH](sA, ss) = y(oo,cx}) . (A.4) 

aA-»0+ Jo- Jo- SB-.CX) 



We then also define [cf. Eqs. 



/ PA(ys)dys=£[p,](oo,0+) = y(oo,0+), (A.5) 
Jo 

i- oo 

/ Ps(yA)dyA = £[p,](0+,«)) = Y(0+,«)). (A.6) 
Jo 



Pa 
Pb 

Using Eq. (|45[), we find Pa = exp(— ptta), Pb = exp{—pnB), and Pab = exp(— pttaus), where tta is the area of the support of 
Wa, ttb is the area of the support of wb, and ttaus is the area of the union of the two supports. This result is of course not 
surprising and has been already derived in the paragraph before Eq. (|6^) using a different approach. 

For vanishing weights, we decided to use the following prescription: We discard, in the ensemble average for Cov(/; Ga, 0b) , 
the configurations {On} for which the function / is not defined either at 0a or at 0s. In order to implement this prescription, we 
can explicitly modify the probability distribution py and exclude "by hand" cases where the denominator of Eq. ( p^ vanishes; 
for the purpose, we consider separately cases where wa or wb vanish. We define a new probability distribution for (yA,yB) 
which accounts for vanishing weights: 

\Py{yA,yB) ifwA/O, 
- , s , [PyiyA,yB) - PA{yB)5{yA)]/il - Pa) if wa = 0, wb / , 

Py{yA,yB) = ■{ \ N D ^ N //I D^ -f n ("^-^^ 

[Py[yA,yB) - PB(yA)S{yB)\/{l - Pb) if wa / 0, wb = , 

[Py{yA,yB) ~ PA{yB)S(yA) ~ PB{yA)S{yB) + PABS{yA)5{yB)]f if wa = 0, ws = . 
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Order Mij 

ji+j) j-0 

1 - - - - 

10 - - - 

2 ^520 P'S'ii pSo2 - - 

3 pSso pS'ii pSi2 pSo3 - 

4 ^5*40 + 3p^5'2o pSso + 3p^S'2o5ii ^^22 + p^5'o25'2o + 2/9^S'iiS'ii pS'oa + 3p^ 5*02 5*1 1 pSo4 + 3p^So2 



Table B.l. Moments Mij up to the fourth order. The table shows, for each row, the values of Mij with {i + j), the 
order, fixed. Hence, for example, the row for order 2 shows M20, Afn, and Mq2 in sequence. 



We recall that v = 1/(1 — Pa — Pb + Pab)- In constructing this probability, first we have explicitly removed the degenerate 
situations, then we have renormalized the resulting probability. Note that the normalization factor in the last case, namely 
1 — Pa — Pb + Pab, comes from the so-called "inclusion-exclusion principle" (1 — Pa — Pb + Pab is the probability that both 
fA and fB are defined). Using this new probability distribution in the definition (^2|) for Y we obtain 

if lUA / 0, Wfl / , 

Y{sa,sb) = <'„o(^^,^. _,o(s.. oo)i r, ^ ..... . (A-8) 




if IDA 7^ 0, Wb ~ 

Finally, we need to change the normalization factor in Eq. (^) in order to account for cases where yA or yB are vanishing. 
Indeed, the correcting factor C(wa,wb) has been obtained by assuming that all objects can populate all the plane with uniform 
probability distribution [cf. Eq. (|25[)]; now, however, a fraction (Pa + Pb — Pab) of configurations have been excluded. Hence 
we have 

C{WA, wb) = ^_p^J'1^^p^^ ^[Y]{wa,wb) . (A.9) 
This complete the discussion of vanishing weights. 



Appendix B: Moments expansion 

In Sect. ^ we have written the moments expansion for C(wa,ws). Here we complete the discussion by providing a proof for 
that result. 

At high densities, yA and j/s are basically Gaussian random variables with average values yA and i/s (we anticipate here 
that these averages are given by the density p). Hence, we can expand them in the definition of C{wa,wb)'- 



C{wa,wb) — p 



AyA / dys 
Jo 

) r 00 

AyA I AyB 



PyiyA,yB) 



(wA + yA){yJB + yB) 

Py{yA,yB) 



i,j=0 



{wA + yA)(wB + yB) 



E 



VA ~ yA 

WA + yA 



E 

j=0 



VB - yB 
WB + yB 



{yA+WA)^+^{yB +Wb)^ + ^ ' 

where Mij are the "centered" moments of py-. 

r 00 r 00 

Mij = / AyA / AyBPy[yA,yB){yA-yAy{yE 
Jo Jo 



Vb) 



The centered moments can be expressed in terms of the "un-centered" ones, defined as 



AyA / AyBPy{yA,yB)yWB = {-!)'+' Y^^''\o,0) 



(B.l) 



(B.2) 



(B.3) 



Here Y^'"'-'^ (0, 0) is the i-th partial derivative on sa and j-th partial derivative on sb of Y{sa, sb), evaluated at (0, 0). These, in 
turn, can be expressed as derivatives of Q. For the first terms we have 



y<'''°'(0,0) =F(0,0) = 1 , 

y (^'O' (0, 0) =pQ(^'°' (0, 0) + [pQ'^'"' (0, 0)] ' , 

y (1'^' (0, 0) =pO<^'^' (0, 0) + [pq(^-") (0, o)] [pq'-"-''' (o, o)] 

y f"'^) (0, 0) =pQ('''^) (0, 0) + [pQ(°-^' (0, 0)] ' . 



y (0, 0) ^pQ^'^"' (0, 0) , y (o, o) =p(3^°''' (o, o) 



(B.4) 
(B.5) 
(B.6) 
(B.7) 
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Finally, the derivatives of Q can be evaluated as 



(B.i 



where Sij, we recall, is given by Eq. (^). Note that 5oi = Sio = 1 because of the normalization of wa and wb, and thus, as 
already anticipated, i/a = ys = P- In summary, we find 



Moo = 1 , 

Mil = Mil - MioMoi = pSii 



Mio = Moi = , 

Mo2 = M20 - {Moif = pS2 



M20 = M2 



{Miof = pS2 



(B.9) 
(B.IO) 



We stress that, in general, it is not true that Alij = pSij (more complex expressions are encountered for higher order terms; cf. 
the last term in Eq. (pq)). Finally, we can write the expansion of C{wa,wb)- 



C{wa,wb) 



+ 



P^S20 



+ 



p'Sii 



+ 



P^So2 



{p + wa){p + Wb) {p + wa)^{p + wb) {p + WA)^ip + WB)^ {p + WA)ip + Wb)^ 



(B.ll) 



This is precisely Eg. (pq). Using the same technique and a littl e rn ore perseverance, we can also obtain higher order terr ns. I n 
particular. Table B.I reports the moments Mij defined in Eq. (B.2) up to the forth order. This table, together with Eq. (B.I), 
can be used to write an accurate moment expansion of C{wa,WB)- 



Appendix C: Varying weights 

In Paper I we have considered a modified version of the estimator (0) which allows for the use of supplementary weights. 
Suppose that we measure a given field fiO) at some positions 0„ of the sky. Suppose also that we use a weight u„ for each 
object observed, so that we replace Eq. (|l|) with 

f((i\ - E^^l fnUnW{0 - On) 

j2„=iUriw(e - e„) 

For example, if we have at our disposal some err or e stimate (t„ for each object, we might use the weighting scheme u„ — 
in order to minimize the noise of the estimator (CI). 

A statistical study of the expectation value of this estimator has already been carried out in Paper I. Here we proceed further 
and study its covariance under the same assumptions as the ones used for the study of Eq. (0) in the main text. However, since 
one of the main reasons to use weights is some knowledge on the variance of each object, we use a generalized form of Eq. (^): 

{[fn - fiOu)] Ifm. - /(0,n)]) = a\un)S„„, . (C.2) 

Note, in particular, that the variance is assumed to depend on u„ (or, equivalently, the weight is assumed to depend on the 
variance). Similarly to Paper I, we also assume that, for each object n, the weight u„ is independent of the position 0„ and of 
the measured signal /„ , and that each u„ follows a known probability distribution pu ■ 

In Paper I we have shown that the average value of f{0) can be calculated using the equations 



)= d^9 d^p4u)[e-''^"™^(^' - 1] = / p4u)Qx{usx)du, (C.3) 
Jn Jo Jo 

Yx{sx)= eyip[pRx{sx)] , (C.4) 

Bx{vx) = pC[Yx]{vx) , (C.5) 



Rx{sx 



{fx)= d'e dupu{u)f{e)Bx{uwx{e))uwx{0) ■ (C.6) 
in Jo 

We now evaluate the covariance of the estimator (CT) using a technique similar to the one used in Sect. ^ We have 

KJaJb) ^ ^ d 01 duipu{ui)--- d 6n duNPuiuN) '"^ r '} . ■ (C.7) 

' ^" Jo Jo Jn Jo En^"™A(0n) Em^™^s(6',„ 



As usual we consider separately the cases n = m and n ^ m, thus obtaining the two terms Ti and T; 



2- 



T ^ [ A^a r A I \ [ A^a r A ( , [f{Sl) + a\ui)]ulwA{ei)wB(ei) 

^1 ^ 71^ / '^uipuiui)--- d On duNPu{v.N)-r^ -|r -— , (C.8) 

^ Jn Jo Jn Jo iL„UnWA[On)\[l^^UmWB{0m.)\ 

NjN-l) f ^2, . , /",2, r, , , f{0l)uiWA(ei)fie2)u2WB{e2) 

T2 = jj:^ / ^1 / duipu{ui)--- d On duN Pu{uN)-r^ -, . — -y . (C.9) 

A" Jn Jo Jn Jo LEn"il"'4(0n)J [Em""i™-B(0m)J 
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Let us introduce new variables vxn ~ UnWx{On) (with X G {A, B}) for the combination of weights, and let us define similarly 
to Eq. ( p6| ) yx = X]^^=2 Then the probability distributions for vxn and yx can be evaluated using the set of equations 

Pv(va,Vb) = \ \ dupu{u)5(vA-UWA{0))5{vB ~UWB{d)) , (C.IO) 

^ Jn Jo 

poo poc poo poo 
Py{yA,yB) = / dVA2 / dvB2Pv(vA2,VB2)- ■ ■ / dVAN / dVB N Pv{vAN , Vb N)S{yA - VA2 Van) 

Jo Jo Jo Jo 

X S{yB - Vb2 - ■ ■ ■ - vbn) ■ (C-11) 

Again, it is convenient to consider the Laplace transforms of these two probability distributions: 

1 p oo p poo 

V{sa,sb)=£.[pv]{sa,sb) = ^ dupuH / d'ee-''-*"'"^'^)-''«""'«(^' = / p{u)W{usa,usb) du , (C.12) 

^ Jo Jn Jo 

Y{sA,SB)=C[py]{sA,SB) = [V{sa,sb)]''~' ■ (C.13) 
In the continuous limit we define instead 



R{SA, sb) 



p roc poo 

d^e dup„H[e-=-*""-*(^'^=«""«(^' -1] = / pu{u)Q{usA,usB)du, (C.14) 
Jsi Jo Jo 

y(sA,Sfl)= exp[pJ?(sA, ss)] • (C.15) 
Finally, the equivalent of the correcting factor C{wa,wb) [cf. Eq. is, in our case, the quantity 

B{va,vb)^P^ r dyA HdyB , P^^^^'^^^ = pC[Y]{va,vb) . (C.16) 
Jo Jo (VA+yA)(VB +yB) 

The quantity B can be used to evaluate Ti: In fact, we have 

= ^ / d Oi duipu{ui) / dyA / dyB Py{yA,yB)-r r^— rr -—— Y 

A Jsi Jo Jo Jo [u-iWa{0i) + yA\[UlWB(0l) +yB\ 

= - f d^Oi f duiPuiui)[f{0i)+a''{ui)]ulwAi6l)wB{0l)B{uiWA{0l),UiWB{0l)). (C.17) 

pJn Jo 
Similarly, for T2 we obtain 

p p 00 p p 00 

T2 = / d^6»i / duip^Ui) / d^e2 / dU2Pu(u2).fi0l)fi02)uiWAi0l)u2WBi02) 

Jn Jo Jn Jo 

X B{uiWa{0i) +U2Wa{02),UiWb{0i) +U2Wb{02)) ■ (C.18) 

The final evaluation of Cov(/) then proceeds similarly to what done in the main text for the estimator (jl|). 

Appendix D: Properties of the Laplace transform 

For the convenience of the reader, we summarize in this appendix some useful proper ties of the Laplace transform. Proofs of 



the results stated here can be found in any advanced analysis book (e.g. Arfken 1985). Although in this paper we have been 
dealing mainly with Laplace transforms of two-argument functions, we write the properties below for the case of a function of 
a single argument for two main reasons: (i) The generalization to functions of several arguments is in most cases trivial; (ii) 
Several properties can be better understood in the simpler case considered here. 

Suppose that a function f{x) of a real argument x is given. Its Laplace transform is defined as 

poo poo 

£[f]{s)= hm / dx' f{x')e-'^' = dx/(:r)e-^^. (D.l) 



Note that we use 0^ as lower integration limit in this definition. 

The Laplace transform is a linear operator; hence, if ct and 13 are two real numbers and g{x) is a function of real argument 
X, we have C[af + Pg] = aC[f] + l3C[g\. 

The Laplace transform of the derivative of / can be expressed in terms of the Laplace transform of /. In particular, we have 

£[/'](«) = .£[/](.) -/(0-) . (D.2) 
This equation can be generalized to higher order derivatives. Calling the n-th derivative of /, we have 

n-l 

(s) = s"£[/](s) - ^"~"'f^"\0~) ■ (D.3) 
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Surprisingly, this equation holds if, for n negative, we consider to be the — n-th integral of /; note that in this case the 
summation disappears. Hence, for example, we have 



f{x') dx' 



(s) = -mis) 



(D.4) 



Often, properties of the Laplace transform co me in pairs: For every property there is a similar one where the role of / and 
£,[f] are swapped. Here is the "dual" of property (D.2): 



C[xfix)]{s) 



ds 



or, more generally. 



C[xy{x)](s) = {-1) 



ds" 



(D.5) 



(D.6) 



A similar equation holds for "negative" derivatives, i.e. integrals of the Laplace transform. In this case, however, it is convenient 
to change the integration limits to (s, (x>). In summary, we can write 



£[.f{x)/x]{s)^ af]is')ds'. 

J s 

Given a positive number a, the Laplace transform of the function / shifted by a is given by 
C[fix-a)R{x)]{s)^/:[f]{s)e-^^ , 

where H is the Heaviside function defined in Eq. (po|). A dual of this property can also be written: 
C[f{x)e''-]{s)=C[f]{s~b). 

Finally, we consider two useful relationships between limiting values of / and £[/]: 



(D.7) 

(D.8) 
(D.9) 



lim f{x) = lim s£.[f]{s) , 

x^0+ s~>oo 



lim fix) = lim s£[/](s) 



(D.IO) 
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